Dynamic localization in quantum dots: analytical theory 
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We analyze the response of a complex quantum-mechanical system (e.g., a quantum dot) to 
a time- dependent perturbation <f>(t). Assuming the dot energy spectrum and the perturbation 
to be described by the Gaussian Orthogonal Ensemble of random matrices we find the quantum 
corrections to the energy absorption rate as a function of dephasing time t v . If <j>(t) is a sum of 
d harmonics with incommensurate frequencies, the quantum corrections behave similarly to those 
of conductivity 5ad{t v ) for the d-dimensional Anderson model of the orthogonal symmetry class. 
For periodic perturbations, the leading quantum corrections are generically absent as in the systems 
of the unitary symmetry class. Exceptions are the harmonic perturbation and all other periodic 
perturbations <f>(t) that obey the generalized time-reversal condition <j>(—t + t) — <f>(t + r), where r 
is a certain shift. Such cases fall into the quasi- Id orthogonal universality class. 

PACS numbers: 73.21. La, 73.23.-b, 73.20.Fz, 78.67.Hc 



The process of energy absorption by a quantum system 
with a time-dependent Hamiltonian underlies a large part 
of modern physics, both fundamental and applied. A 
generic Hamiltonian can be written in the form: 



H{t) = H + V</>(t), 



(1) 



where we explicitly separated the time-independent 
part Hq and the external perturbation V with the time 
dependence specified by a given function 4>(t). Most of- 
ten the relevant case is that of the classical ohmic Joule 
absorption. The simplest nonlinear effects that restrict 
the absorption rate are the saturation effects that origi- 
nate from an upper bound on the spectrum of H (as in 
the textbook example of a two- level system). 

In the past two decades attention of the scientific com- 
munity was drawn to a different and much less trivial 
example of saturation when the spectrum of Hq is es- 
sentially unlimited in the energy space, yet, after a cer- 
tain time the absorption stops. This so-called dynamic 
localization (DL) in the energy space was observed in nu- 
merical simulations on the kicked quantum rotor (KQR) 
- particle on a circle with Hq = —d^JdO 2 and <p{t) be- 
ing a periodic sequence of <5-pulses [|lj, as well as in an 
actual experimental realization of the KQR - trapped 
ultracold atoms in the field of a modulated laser stand- 
ing wave 0]. The mapping of the KQR to the quasi- 
random Id Anderson model has been done in Ref. ||, a 
similar analogy was exploited by Gefen and Thouless Q 
to demonstrate the DL in a mesoscopic disordered ring 
threaded by a magnetic flux growing linearly in time. In 
Ref. pi an analogy between the KQR and band random 
matrices was pointed out, the latter have been reduced 
to a Id nonlinear a model [|| . In Ref. ffl the direct corre- 
spondence between the KQR and a Id nonlinear a model 
was demonstrated. 

On the other hand, numerical simulations for a 5- 
kicked particle in an infinite potential well that differs 



from the KQR only by the boundary conditions, revealed 
no DL (§). This example shows that DL is not a conse- 
quence of the one-dimensional character of the energy 
space but it depends on the details of both the unper- 
turbed system and the perturbation. It is clear that the 
results on DL || obtained in the framework of the ortho- 
dox kicked rotor (standard map) model jlj] cannot be au- 
tomatically applied to a generic quantum mechanical sys- 
tem under arbitrary time-dependent perturbation. Even 
less clear is the status of a peculiar KQR model with the 
time-dependent perturbation characterized by three in- 
commensurate frequencies where the phenomenon similar 
to the 3d Anderson localization-delocalization transition 
has been found numerically 

The most general assumption about a complex 
quantum-mechanical system would be the randomness 
of the Hamiltonian. It is well-known M] that at energies 
smaller than the Thouless energy E c description of com- 
plex quantum systems falls into one of the three univer- 
sality classes, each corresponding to a Gaussian ensemble 
of random matrices. Considering a non-magnetic elec- 
tron system with the spin-rotation symmetry, we arrive 
at Eq. (|l|) with H and V from the Gaussian orthogonal 
ensemble (GOE). An analogous Hamiltonian was consid- 
ered in connection with laser pumping of complex organic 
molecules [[ll], but DL was not a target at that stage. 
The problem of DL in systems described by the Hamilto- 
nian (|l|) has been addressed in Ref. Jl2fl . However, most 
of results of this study are qualitative in character and 
DL has been demonstrated numerically only for the case 
where the typical amplitude of harmonic perturbation 
V and its frequency to are of the order of the mean level 
spacing 8 for the unperturbed system. This numerics has 
also shown that DL is unstable with respect to adding 
even a small amount of noise in the time dependence. 

In this Letter we develop an analytical theory of DL 
for a closed system described by the Hamiltonian (fy) 
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with Hq and V from the GOE of random matrix theory 
(RMT) and a general time dependence <j>{t). We show 
that in this model the existence of DL and even the clas- 
sical ohmic Joule regime with constant absorption rate 
depends strongly on the character of <fi(t). For periodic 



^t) = En Ant 



with \A n \ decreasing as Ifn 6 or 



slower, the absorption rate diverges and the ohmic regime 
does not exist in the limit of infinite matrix size N. In 
particular, this happens if <f>{t) is a periodic ^-function as 
in the orthodox KQR model. This result alone demon- 
strates the difference between the KQR and the time- 
dependent random matrix theory. 

For a periodic <fr(t) with rapidly decreasing A n , we ob- 
tain that the ohmic Joule absorption rate Wq is modified 
by the quantum interference effects responsible for the 
DL. We calculate the first (one-loop) quantum correc- 
tion, which determines the weak DL, and show that it is 
only nonzero if the condition 



-t + r) = (f>{t + T) 



(2) 



is fulfilled for a certain choice of r. In this case the correc- 
tion 8W(t) oc \ft grows with the time t of the action of the 
time-dependent perturbation, like in the quasi-ld Ander- 
son localization in the orthogonal symmetry class. At a 
certain time scale t* the first correction 6W(t*) becomes 
comparable with the Joule absorption rate Wo indicat- 
ing a crossover from the weak to strong localization in the 
energy space. On the other hand, if dephasing time t v 
is shorter than t», then dissipation remains ohmic, with 
the rate being smaller than Wq by the value of SW(t v ). 

The condition (|^) is a generalization of the time- 
reversal symmetry condition cj>(t) — <fi(—t) for the case 
when a shift of the time origin does not matter jOI . For 
periodic functions with several harmonics Eq. (0) can be 
fulfilled only for special choice of relative phases. How- 
ever any pure harmonic function obeys this condition and 
thus the monochromatic perturbation is rather an excep- 
tion than a paradigm of a periodic perturbation. For 
generic periodic functions which do not obey Eq. (||) the 
quantum correction emerges in the two-loop approxima- 
tion, as in systems of the unitary symmetry class. 

Finally, we find that for <j>(t) being a sum of d har- 
monic functions with incommensurate frequencies the 
first quantum correction to the absorption rate is sim- 
ilar to the weak localization correction to conductivity 
of a cZ-dimensional disordered system of the orthogonal 
symmetry class. This gives an analytical support of the 
numerical observation of Ref. |B for the KQR model. 

Qualitative picture and mapping. — The similarity and 
difference between DL and AL, as well as between KQR 
and RMT can be qualitatively understood using the ex- 
act correspondence between a quantum system under 
a (multi-)periodic time-dependent perturbation and a 
tight-binding lattice model with the time-independent 
Hamiltonian. Consider a system with energy levels ("or- 
bitals") Ei under a harmonic perturbation Vwe~ luJt + 
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FIG. 1: Lattice analogy for a monochromatic perturbation. 



Vj*;e iw *. As follows directly from the Schrodinger equa- 
tion, its time evolution can be expressed in terms of eigen- 
functions and eigenvalues of another system (see Fig. [l]), 
obtained from the original one by replicating it into a 
one-dimensional lattice and shifting the levels of each 
consequent site by fao so that the energy of the "or- 
bital" I on the site s is given by Ei >s — E\ — sfaj, and 
introducing the coupling between the Ith orbital of the 
site s and the I'th orbital of the site s + 1 by the station- 
ary perturbation matrix element Vw . In the same way 
one can show that higher harmonics in the perturbation 
(2oj, 3oj . . .) correspond to the coupling to next neighbor- 
ing sites (second, third, etc.), while the presence of sev- 
eral incommensurate frequencies u}\,...,u}$ requires a d- 
dimensional lattice with sites s = [s\, .. . Sd), on-site en- 
ergies Ei — S\UJi — ... — s^Wd, and the matrix element 
at the ith frequency corresponding to the coupling along 
the ith dimension. 

The d-dimensional space with "orbitals" on different 
"sites" coupled by the matrix elements of perturbation, 
arises quite naturally from the notion of quasi- energy p[ . 
However, this model is far from the conventional Ander- 
son model, where there is only one orbital per site with a 
random energy and from the Wegner's ./V-orbital model 
Jl5| where there is a fixed number of orbitals with random 
energies and random hopping integrals. The point is that 
being a copy of the original spectrum, the set of orbitals' 
energies are not independent at different sites. Another 
difference is that the number of orbitals is infinite in the 
relevant limit of the infinite matrix size. 

From the above mapping it is clear why slowly decreas- 
ing amplitudes of harmonics kill localization in RMT. 
It happens because of the possibility of long-range hops 
over the sites. It is also clear why the long-range hops 
are not so dangerous for localization in the KQR case. 
The point is that in the orthodox KQR model, the time- 
dependent perturbation is coupled to cos 9 which in the 
basis of eigenfunctions e lW of the unperturbed system 
corresponds to matrix elements Vid±i that may connect 
only neighboring orbitals. For long enough distance be- 
tween sites the coupled orbitals run out of resonance and 
the hopping is suppressed. Should the basis change to 
sm(W) as in the case of infinite potential well in Ref. [0, 
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the matrix elements Vi t i' oc l/\l — l'\ become long-range 
in the orbitals' space, so that the resonance hopping be- 
tween remote sites may still occur, though only between 
remote orbitals. For the RMT case Vn' does not de- 
crease with the distance |Z — l'\ and the condition for the 
long-range resonance hopping is further improved. 

This consideration shows that the KQR with <5-kicks 
and the time-dependent RMT with harmonic <j>{t) are 
in fact two complementary models. In the former re- 
mote sites are connected but only for neighboring orbitals, 
while in the latter remote orbitals are connected only on 
neighboring sites. Sometimes - as in the case of sev- 
eral harmonics with incommensurate frequencies - this 
complementarity makes two models equivalent but this 
equivalence is fragile and needs to be checked in every 
particular case. 

Description of the formalism. — Consider a closed sys- 
tem of non-interacting fermions with the single-particle 
Hamiltonian (|l]), Hq and V being real symmetric random 
N x N matrices distributed with probabilities 



Vh oc exp 



7T 2 trff 2 
4NP 



Vv oc exp 



wtiV 2 



4T5 



respectively, 6 being the mean level spacing at the center 
of the band, and T measuring the sensitivity of the energy 
levels Ei to the variation of <j>: ({dEi/ dcf) 2 ) = (2/n)T8. 
Such a form of time-dependent random matrix theory 
corresponds to a quantum dot under a perturbation with 
characteristic frequencies lu <C E c p6[ . 

To treat the dynamical problem specified by the Hamil- 
tonian (|l|), we employ the nonlinear a- model approach re- 
cently developed in Ref. fl^| o n the basis of the Keldysh 
nonequilibrium formalism |l8|, [l9| . In the limit N — > oo, 
the effective action of the er-model reads: 



S[Q] 



26 



tv{iT 3 a Stt'd t Qtt'} 



dtdt'. (3) 



Here Qu> is a 4 x 4 matrix in the direct product of the 
2x2 particle-hole and Keldysh spaces. Pauli matrices in 
these spaces are denoted by Tj and <7j, respectively. The 
first term in S is the standard random-matrix action j2(J 
responsible for the whole spectral statistics, while the 
second - kinetic - term accounts for the effects of the 
time-dependent perturbation. The matrix Q is subject to 
the constraints (Q 2 )te = TWo^t*' and T20\Q<J\T2 = Q T , 
where the product and the transpose involve the time 
arguments too. The saddle point of the action (||) is 



A, 



(o) 



Stt> 2F£ 







73, 



(4) 



where the function F^) satisfies the kinetic equation 

[d t + d t/ + rm)- ( j ) (t')) 2 ]F t { ° ) = 0. (5) 



F u > is the electron distribution function in the time rep- 
resentation; a more familiar quantity is its Wigner trans- 
form FE(t) = J dre lET F t+T /a t— t /2- I n equilibrium with 
temperature T it is F E = t&nh(E/2T). Out of equilib- 
rium it satisfies the Wigner-transformed Eq. (|B|), which 
after averaging over fast oscillations in t reduces to the 
diffusion equation in the energy space: 

p(o) 



[d t - d d%] Fg> (t) = o, d = r (d4>/dt) 2 , (6) 

the over line meaning the average over time. Equation (^) 
gives ohmic Joule absorption rate 

W a = D/S. (7) 

The saddle-point expression ^ is valid provided that 
(i) the perturbation is sufficiently fast, d t 4> S> 5 3 / 2 r _1//2 , 
which is the anti-adiabaticity condition ensuring that the 
spectrum is quasi-continuous and (ii) interference effects 
responsible for DL are neglected. 

The perturbative correction to Eq. (Q) for the case of 
linear bias <fi{t) = vt was calculated in Ref. [j!?]]. There, 
the interference is ineffective and quantum corrections 
describe the high-ii tail of the crossover from the Kubo 
to the Landau-Zener regimes of dissipation Con- 
trary, in this Letter we concentrate on the case of large 
d t <f> > 5 3 / 2 r -1 / 2 but bounded <j>(t), when the saddle- 
point approximation is invalidated by interference effects. 

Quantum corrections to the mean-field absorption 
rate (0) can be obtained in the regular way by expanding 
over Gaussian fluctuations near the saddle point. The de- 
viation of the Q matrix from A is parametrized Jl7| by the 
diffuson and cooperon modes btv and aw with the bare 
propagators: (6 t+f _ 6^, t , ) = {25/tt) 5{r) — rf) ^(t, t') and 
(a* + t„a t * f ) = {6/-!r)6(t-t')C t (r),r]'), where 



V n (t,t') = 6(t-t') exp 
Ct(v,v') = 9(v ~ V') ex P 



%i(r)dT 

t' 

1 rn 

2 J v > 



(8) 



(9) 



and we have denoted t± = t ± rj/2, t'± = t' ± rf /2, and 

7l? W = r[0(t + )-0(t_)] 2 . 

In the presence of fluctuations, the average matrix (Q) 
still has the form (|]) but with the saddle-point F^ 
substituted by the renormalized electron distribution F 
which determines the energy absorption rate: 

W{t) = d t (E(t)) =-j Kmd t d n F t+v/2!t _ v/2 . (10) 

The one-loop quantum correction to (Q) contains a 
diffuson and one cooperon loop coupled by a Hikami box. 
We evaluate this diagram for a generic perturbation <f>(t) 
switched on at t = 0. Using the asymptotics F t+ t_ ~ 
(iirr)) -1 following from the property liruE_>±oQ Fe {t) = 



±1, we obtain for the absorption rate: 



r 



W(t) = W + - / 0W(t-OC t _ 6/2 (£,-Ode- ( n ) 



The second term in Eq. ( |11[ ) is the one-loop quantum 
interference correction to the ohmic Joule heating (Q). 
Expression ( )ll| ) is the main result of this part of the paper 
and will be the base for the subsequent considerations. It 
can be also obtained from the conventional diagrammatic 
technique [|l6| p2[ . 

Results. — First, we consider a periodic perturbation: 
<t>(t) = ^2 n A n cos(ncut—ip n ), with the diffusion coefficient 
D = (lu 2 T/2) J2 n n 2 A 2 n . If A n decrease as rr z l 2 or slower 
(e. g. for 5- kicks), the diffusion coefficient and dissipation 
rate diverge. This divergency is cut either by a finite 
size N of the matrix or by the interaction effects that is 
beyond our analysis. 

To study the long-time, period-averaged dynamics at 
t, £ 3> 1/ui we can approximate 



C t -£/ 2 (£,-£) ~ e 



- 2r i E„ A ™ sin 2 [raw(t-{/2)- V „] 



(12) 



For a particular choice of phases, ip n — nip, there exists a 
set of points = 2t — 2{<p + nk)/uj with integer k where 
the cooperon ( |I2| ) is equal to unity (Fig. ||a). Existence of 
such no-dephasing points ]l3| is equivalent to the general- 
ized time-reversal symmetry (Q) of the perturbation. At 
large £ only £ m & contribute to the integral (11) (other- 
wise the cooperon is exponentially small) which then can 
be calculated with the steepest descent method. Per- 
forming summation over the no-dephasing points we 
obtain a negative and growing in time quantum interfer- 
ence correction to the ohmic absorption rate ([?]): 



W(t) 
Wo 



= 1 - 



2S 2 



(13) 



where n 2 = J2 n 71,2 A n and the limit t >■ 1/uj, 1/T is im- 
plied. The \ft dependence is remarkably similar to the 
y/T^ dependence of the first quantum correction to the 
conductivity of a particle in a quasi-one-dimensional dis- 
ordered sample with the phase relaxation time t v f23]| . In 
our case the relative quantum correction becomes com- 
parable to unity at time i*. This is an indication of DL 
at the characteristic energy scale E* ~ y/Dt* ~ u>Tn 2 /d. 
Taking 6 = 3 ^eV, Hu = 40 ^eV (lu/2it w 10 GHz), 
Tn 2 = 10 neV (corresponding to the microwave electric 
field of a few V/m over the dot size of 1 fim) we obtain 
U ~ 10 ns, E* ~ 400 /xeV ~ 5 K. 

When the generalized time-reversal symmetry (0) is 
absent for <fi(t), the dephasing rate is positive and sep- 
arated from zero by a finite gap (Fig. ||b). In this case 
the integral over £ converges exponentially, so the first 
quantum correction stays small even at t — > oo, as in the 
systems of the unitary symmetry class. 

Next, we consider the case of d incommensurate fre- 
quencies ui n ; 4>{t) = Y^ n= iA n cos{uj n t — ip n ), for which 




FIG. 2: Perturbation-induced dephasing rate [the coefficient 
at £ in the argument of the exponential in Eq. ([[i])] vs. £: 

(a) periodic 4>(t) obeying Eq. (g): a regular array of zeros; 

(b) general periodic 4>(t): a gap; (c) quasi-periodic </>(i) with 
two incommensurate frequencies: neither a gap nor a regular 
array of zeros. 



D = (r/2) J2 n =i w n^n- For incommensurate frequencies 
the relationship between the phases ip n does not matter. 
The reason is that in this case there can be only one 
no-dephasing point at most which cannot make a deci- 
sive contribution to the integral, in contrast to the pe- 
riodic case with an infinite number of such points. On 
the other hand, in the incommensurate case the dephas- 
ing rate can become arbitrarily small (pseudo-gap) in an 
infinite number of points regardless of phases (Fig. ||c), 
and the integral is dominated by these points. 

To calculate the first quantum correction in Eq. ([ll]) we 
use the cooperon ([l2]) with nuj substituted by ui n , expand 
it into a d-dimensional sum involving the modified Bessel 
functions I v (z), and average over t which is significantly 
simplified due to incommensurability. The analytical ex- 
pression for the case of arbitrary relation between the 
amplitudes A n is bulky. We write the resulting expres- 
sion valid at t 3> 1/w, 1/T for the case when all A n = 1: 



W(t) 
Wo 



s 

7rr 



dz e 



- z d 



[*>(*)] 



d-1 dI o( z ) 
dz 



(14) 



Using the asymptotic form Io(z) ~ dIo(z)/dz ~ e z j \/2~nz 
at z » 1, we recover Eq. (13) for a monochromatic 
perturbation. For d = 2, the relative correction is 
(5/2tt 2 T) lnTi, whereas for d > 2 it saturates oc t 1_d / 2 in 
the limit of large t, in complete analogy with the behavior 
of the quantum correction of Ref. |23| in d dimensions. 

In conclusion, we have developed an analytical ap- 
proach based on the zero-dimensional, time- dependent 
nonlinear sigma-model and obtained the weak dynami- 
cal localization in complex quantum systems under time- 
dependent perturbation described by the random matrix 
theory. The character of energy absorption in such sys- 
tems is determined entirely by the frequency spectrum 
of a time-dependent perturbation. In particular we ob- 
tained no DL for the time-periodic (5-function perturba- 
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tion, and the dynamical localization corections similar to 
the d-dimensional weak localization corrections to con- 
ductivity if the perturbation is a sum of d incommensu- 
rate harmonic functions. 
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